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Abstract 

We study a homogeneous linear second-order difference equation with 
nonconstant and noncommuting operator coefficients in a vector space. We 
build its exact resolutive formula consisting of the explicit noniterative ex- 
pression of a generic term of the unknown sequence of vectors. Some non- 
trivial applications are reported in order to show the usefulness and the broad 
applicability of the result. 



1 Introduction 

Difference equations naturally occur and play a central role whenever the prob- 
lem under scrutiny or the phenomenon under investigation allows a mathematical 
formulation traceable back to the set of natural numbers. A broad variety of sit- 
uations in biology, economics, dynamical systems, electrical circuit analysis and 
other fields |1, chap 1] is indeed modeled by difference equations often demand- 
ing, however, quite different resolution methods due to peculiar aspects of their 
mathematical nature. Thus, like in the case of the differential equations, a useful 
classification of the difference equations has emerged too. 

The classical field of difference equations deals with linear or nonlinear equations 
where the unknown is a real or complex-valued function defined on a countable 
domain which in turn may be thought of as coincident with N* without loss of 
generality. These equations are conveniently classified as scalar difference equa- 
tions [:2, chap 2,3,4]-[|6l to discriminate them from the so-called matrix difference 
equations 0, ||8l where the unknowns constitute a sequence of d x d matrices with 
entries in C as well as possibly all the given coefficients. 

Quite recently a new class of difference equations has been introduced, namely 
the class of the linear operator difference equations ||9]- The peculiarity of such 
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equations with respect to the matrix difference equations is that in this case the 
unknowns are elements of a given abstract vector space V whereas the "coefficients 
"are linear operators acting on V. We clarify this important point by introducing 
the following linear second order operator difference equation we are going to 
investigate and solve in this paper: 

Y n+2 = £ (n)Y n + £i(n)Y B+ i, n € N*. (1) 

This equation lives in a vector space V over a field F and the nature will not be 
further specified. The null vector of V is denoted by 0, while £o(n) and L\ (n) , n S 
N are two families of linear operators acting in V. To keep the investigation at a 
general level as much as possible, we do not introduce any constraint concerning 
the commutability of these operators in and between the two classes. 

The aim of this paper is to build the general solution of equation £T|) when its 
coefficients are nonconstant and noncommutative. This means to give a nonitera- 
tive expression of Y n , for a generic n, in terms of an appropriate ordered sequence 
of operators Co(n) and C\(n). Of course this expression will contain two summa- 
tion constants to be specified starting from given initial conditions. 

The main result will be derived by mathematical induction in Section 3 and 
some applications will be presented in the successive two sections. 



2 Mathematical preliminaries 

We start by recalling that when Vre € N* £o( n ) = £o an d C\{n) = C\ the notation 



u + v 



possible dis- 



{Cq^C^} with u, v G N* expresses the sum of all the 

tinct permutations of u /^-factors and v £i-factors, where m := mm(u,v) ifTOll . 
For example {C^C^} = CqC\ + C\Cq and we define {C^C^} := I putting 



consistently 






(u) r (v)i 



1. It is possible to convince ourselves that any additive con- 



tribution to {C^'jC-Y'} may be represented as a formal product of r(r > 1) powers 
of £q alternated with r powers of C\, of the form 



rT\ ns\ pT2 rS2 rr r ps r 
Lj \ ^0 • • • ^-O ■ 



(2) 



The parameter r, hereforth called the length of the expression ©, is an integer 
number running from 1 to m + 1 = tm- The 2r integers exponents n, . . . , r r and 
s \ , . . . , s y are numbers fulfilling the following conditions 



i=l 
r 

i=l 

Si > 1, T\ > 0, T r > 1, S r > 



U 



V 



(3) 



n >i,Si> l, 



2,...,r-l 
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Hereafter the symbol v represents the set of all the possible pairs of r-uples 
f = (n, . . . ,r r ) and s = (si, . . . , s r ) satisfying eq. ©. Expanding the symbol 
{C^C^} in accordance to its definition we get contributions possessing all the 
compatible lengths and for any length all the possible terms obtained in accordance 
with eq. ©. It is possible to show that the number of additive contributions of 

equal length r to {C^£^} is given by f ^ U ^ ^ f ^ V ^ ^ so that using the 

Vandermonde identity ifTTTl 

£(i-l)(m + l- i) = ( m ) (4) 
we derive the following length-based representation of {C^C^} 

m+l 

{4 U) 4^} = ^ J2 (5) 

where the summation 

Co^i £l 2 C s 1 >...£Z£ s {- (6) 

(f,s)<=Sr v 

runs over all the pairs (f, s) compatible with the conditions (f3]). The operator co- 
efficients appealing in eq. (fl]), depend on n and do not commute with each other 
in general. In such a case powers like Cq or become ambiguous as soon as t\ 
or sj exceeds 1. Thus, to keep as much as possible the advantages of the notation 
adopted in expression ©, we introduce the following descending products of op- 
erators 

Pr 

]J X p := X( Pr ) . . . . . . X( Pl ), Pl G N*. (7) 



p=pi 



Pr 



postulating that [ X p = I when p r < p\. For the sake of clarity, we stress 

p=pi 

that the symbol introduced in the left-hand side of eq. © singles out the inverted 
ordered product of the noncommuting p r operators X pi , . . . , X Pr . 

Exploiting eq. |7]) we get rid of the ambiguity possessed by expression (O 
due to the noncommutativity of d, i € {0, 1} by substituting it with the following 



3 



ordered operator [Cq 1 C^ 1 C t 2 C S2 . . . C^ r C s {] q defined as follows: 

\ri~i rSi fT2 rS2 rr r fSr] . 

l~0 '-I ^0 L -\ ■ ■ ■ ^0 M \q -~ 

2ii 



T \lC [k q 

il=0 



r 2 -l 

]J C ( k q - 2ti - si - 2i 2 



12=0 
7> — 1 



"ll £i( k q -2n - 

Ji=0 

U Ci[k g - 2(n + r 2 ) - Si - J 2 

i2=o 



i r =0 

where 



U £ { ^-2 T £ 1 Ti- S £ 1 Sz -2i r 
z=i z=i 



s r — 1 



k q := 2u + v 



jr = 



9 = 1,2. 



i=i 



(9) 



In the next section we shall see that the peculiar "order-explosion"of eq. (O as 
given by eq. © provides a very useful tool to write down the solution of a generic 
Cauchy problem associated to eq. (fl]). In the following the symbol {C^ C^} q 

represents the sum © where each additive contribution {C U ^ C^} is replaced by 
the associated ordered operator that is 



m+l 



(10) 



To clarify the notation, we develop for example {£q>C^} 



{4 2) 4 1} i 



i- 



EV- \ rT\ ps\ r>T2 rS2 rr r rSrl 
L*-0 L -\- ■■• ^0 Lj l Jl 

r=l (f ,s)e55 a 

Since fci = 4, the three terms acquire the following explicit form: 



(11) 



U Co 

ii=0 



2u 



U A ji 

ii=o 
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o 

U A 

h=o 



Ji 



U £ 

12=0 



(8) 



£ (4)A(2)£i(0), 

1 = A(4)£ (3)£o(l), 



[4£l4£?]i = £o(4)A(2)A)(l). 
We thus conclude that 

{£§£!>! = >Co(4)£ (2)A(0) + A(4)£o(3)£ (l) + 4,(4)A (2)4,(1). 
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3 The resolution of a Cauchy problem associated with eq. 

It is convenient to extend the definition of the operator {Cq £±}q to negative 
integer values of u and v, simply putting, when u or v are negative integers, 
{4" ) 4" ) } 9 ^ = 0,for any Y which Cq, L\ can legitimately act on. We are now 
ready to prove the following theorem constituting the main result of the paper: 

Theorem 1. The solution of the equation ([7]) given that Yq = and Y\ = B may 
be written down as 

Y n = {^r X - 2t }iB. (12) 

t=0 

Proof 

We prove (fl2l by mathematical induction. For n = and n = 1, eq. (fT2l gives 
Yq = and Y\ = B as expected. It is immediate to deduce Y<i = B from both eq. 
® and eq.CE!>. 

Let's suppose now that the first (n + 1) terms of the sequence solution of eq. 
CD) are representable by eq. (IT2l . We must prove that whatever n > is 

[*fi] 

y n+2 = £ {£* ) ^ +1 ^ 2 '} lJ B (13) 
t=o 

r n+l 1 , 

I a I m'+l 

= E E E [^£? 1 ...>cs-/:;i 1 B 

t=o P =i {f,s)eSl n+1 _ 2t 

where m' = min{t, n + l — 2t} satisfies eq. £[]). 
To this end we start by appropriately transforming 

[^\ 

C Q (n)Y n = Coin) £ {CqC^-^B (14) 
t=o 

= C (n) £ {4 _1 ^ +1 - 2t }lS 
t=i 

1-2 J m"+l 

= A(n) EE E [^^...^liB 

t=i r =i (f,s)esr_ 1 , n+1 _ 21 

where m" = min{t— 1, n+l — 2t} < m' and we have used the fact that pi-^]+l = 
pr^]. In the following we are going to show that for eveiy r, 1 < r < m" + 1 the 
following operator relation holds: 

Co(n) £ [CgC? ...££Xrh d5) 

= r E [£?£r..£?xr]i 

(r,s)GS[ n+1 _ 2t/Ti>0 
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where the right-hand side summation symbol means that the ordered operators 
[C&l 1 . . . £%C%\i are not included. 

We point out that Cq(ti)[Cq 1 C^ 1 . . . £Q r £^ r ]i generates permutations of the same 
length r as £q 1 £* 1 ■ ■ ■ C T § C s { with t\ at least 1 necessarily. In addition we stress 
that the range {1,2,..., m" + 1} of r is compatible with the pair u = t,v = 
n + l — 2t. First, using eq. © and taking into account that in this case k\ = n — 2, 
the left-hand side term of eq. (fT5T ) may be developed as follows: 



Co(n) _ _ £ 

(r,s)6SJ'_ 1 n+1 _ 2t 



ri-1 

U C [n-2-2i 1 



ii=0 



si-1 

U A ( n - 2 - 2ti - ji 

51=0 



ll* £ ( n - 2 - 2n - s x - 2i 2 



12=0 



S2-1 

II £i(n-2-2(n+r 2 )-si-i 2 

j 2 =o 



r r — 1 



i r =0 



Sr-l 



II £ n-2-2 £ n - £ *i-2i, 



Z=l 



JV=0 



s r _i 



II £i n-2-2^T| - £ 



(16) 



i=l 



Shifting up the running index i\ by 1 that is putting 1\ =: i\ + 1, we transform the 



first £o-segment 



n-i 

JJ £ (n-2-2i x 

ti=0 



appearing in eq. ([TBI ) into ] [ Co\n 



h=i 



2I\J , which in turn may be rewritten as \J Co ( n — 2I\ J by multiplication with 

^(re- 
introducing the positive index T\ =: T\ + 1 the previous expression ([TBI may 
be cast in the following form 



(?,s)esr 



t,n+l — 2t/Ti=Ti 



Ti-1 

II £ (n-2J 1 

/i=0 



si-l 

II Ci[n-2T 1 - n 

ii=o 



r 2 -l 

JJ £ (n-2T 1 -s 1 -2i 2 



12=0 



S2— 1 



JJ C 1 [n-2(Tx+T2)-sx-j 2 

32=0 



T r — 1 



i r =0 



JJ Coin - 2T X - 2 £ n - £ s r - 2i, 



1=2 



r=l 



S r — 1 



j r = 



(17) 



JJ £i n - 2Ti - 2 £ 77 - £ S[ - j, 



1=2 



where each ordered contribution still has length r and the passage from t — 1 to t 
arises from the "absorbtion"of Co{n). Eq. (fTTl) is exactly the development of the 
right-hand side term of eq. (TT5T ). where T\ plays the role of t\ and then eq. (fT4l) 
may be given the following aspect: 



[^]m"+l 

Co{n)Y n = £ £ 



£ 



i=l r =l (r,s)£S 



^o 1 ^ 1 ...£S r £i r ]iB (18) 



t,n + l-2i/r 1 >0 



To establish a connection between eq. (TT8T ) and eq. (fT3T ) we extract from the latter 
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expression the ordered operators beginning with jCq 1 , with t\ > 



r n+1 1 , 

[-5-] m'+l 

s E E 



(19) 



i=l r =l (f,S)£Sf, 



,+1 — 2t/ri >0 



Here i = is excluded since it implies m' = and then necessarily t\ = 0. If 
m' = min(t, n + 1 — 2t) = t the coiTespondent highest length ru = t + 1 is 
indeed certainly incompatible with the condition t\ > 0. Thus we are legitimated 
to substitute m' with m" = min(t — 1, n + 1 — 2t) in the expression (fT9l . These 
arguments prove that £o(ra)Y^ captures all the ordered operators appearing in eq. 
(fT3T ) effectively beginning with an £cr se g men t (ti > 0). 

In the following in view of eq. £[]) we are going to prove that Ci{n)Y n+ \ coincides 
with the sum of all the other contributions in eq. ([T3T l namely all those ordered 
operators beginning with Cq 1 (n) with n = 0. To this end we exploit the inductive 
hypothesis writing down the action of Ci(n) on Y n+ i as follows: 



d(n)Y n+l = din) EWl~ 2t }iB 
t=o 



(20) 



[f ] m"'+l 

= A(n)E E E [/^.../^ub. 

t=0 r=l (f, S -)e5[ in _ 2t 

where m'" = min{t, n — 2i}. 

We firstly consider the terms of increasing length r in Y n+ \, having n = 0. 
Taking into account that k± = n — 1 when r = 1 we get from eq. d8j the unique 
contribution (t = 0) 



v ' 7 U,ri/r^ —0 



(21) 



£i(n) 



n-l 

IJ A(n-l-ja 

ii=o 



IJ A(n-Ji 

Ji=0 



_j 1= i 
{£g£? +1 }i 



U A U-Ji 



where we have replaced j\ by J\ — 1. This operator applied to is present in eq. 

CES, (t = 0). 

Now we are going to prove that a generic term of Y n+ \ having length r > 1 and 
still with 7~i = satisfies the property 



£i(n) __ £ ...£^1! = 

( T ' S ) e5 t,n-2t/T 1 =0 

E [/^...^rfh. 

( T,S ) e ^,n + l-2t/T 1 =l) 1 s 1 >l 



(22) 
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Expanding the left-hand side of eq. (1221) we indeed get 



A(n)__ £ 

(T,s)G5 t r _ n _ 2t/ri=0 



ri-1 

LJ A ( n - 1 - 2*i 

ii=0 



U A n - 1 - 2ti - ji 

ii=o 



(23) 



T2 — 1 / \ S2 — 1 

]J A I n — 1 — 2t\ — si — 2i 2 

i 2 =0 \ 



]J A ( n - 1 - 2(n + r 2 ) - si - j 2 

J2=0 



r r — 1 / r r _i 

LI A n-l-2 E 71- E aj-2i, 

i r =0 V J=l J=l 



which, putting ji + 1 = J±, s± + 1 = S%, yields 



6> — 1 / T r S r —1 

U A U-1-2Et(- E a« 
j,=o V 2=1 j=i 



E 



-2t/-r 1 =0,s 1 = S 1 >l 



71-1 

]J £ ( n - 2ii 

ti=0 



Si-l 

U A(n-2n-Ji 

Ji=0 



(24) 



T lJ 1 £o(n-2r 1 -5i-2i 2 

12=0 



r m — 1 / Tr-1 Sy-l 

]J £ n - 2 E ti- Si - E si-2i, 

i m =0 \ 1=1 1=2 



S2-1 



U A(n-2(r 1 +r 2 )-5i-j 2 
i2=o 



s r — 1 / r r «i — l 

IJ A n-2£7i-Si - E *J 

jV=0 V 1=1 9=2 



Eq. (1241) is the expansion of the right-hand side of eq. (1221 where the increase from 
(n — 2t) to (n + 1 — 2t) is due to the "absorbtion"of C\{n) as given by eq. (l24l) 
and the condition si > 1 directly stems from the fact that by construction r > 1 
in eq. (1221) . Summing both members eq. (I22l in view of eq. (I20l i over r > 1 and 
t>l(t = 0^r = l) yields 



[f ] m'"+l 
AHE E 



E [£^£^...£^11 



t=l r=2 (f,s)GS. 

v ' 7 t,n — 2t/T^— 

[f ] m"'+l 

EE E [^^...^^ili 

t=l r =2 (r, S -)e5-„ +1 _ 2t/Ti=0si>1 



(25) 



Note that since m!" < m! each ordered operator of length r > 1 compatible 
with {u = t, v = n — 2t} generates an ordered operator compatible with {u = 
t,v = n + 1 — 2t} increasing si by 1 thus rendering in the latter case s\ > 
1. This circumstance means that all the ordered operators appearing in eq. (|25T ) 
determine only the set of all ordered operators having r > 1 and compatible with 
the maximum length m' + 1 under the constraints n = and si > 1. Thus putting 
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together eq. d2Tb and eq. (1251 we claim to have proved that 

[f] m"'+l 

A WE E E [C£C 8 1 1 ...Cl r £i r ]i = (26) 

t=0 r=l (r, S -)6^„_ 2t/Ti=0 

EE E [/^...rjxfli. 

t=0 r=l (r >S -)eS; jn+1 _ 2t/T1=0)S1>1 

where [S] may be substituted with [ 1J -^-] exploiting the fact that for any n odd the 
highest value of t gives rise to contribution for which t\ > 0. We concentrate now 
on the action of £\{n) on the generic terms of Y n+ \ in eq. d2Qb for which t\ > 0, 
that is 



A(n) _ E 

( T ' S ) eS i%-2t/T 1 >0 



(27) 



To this end we observe that under the constraint t\ > the ordered operators as 
given by eq. (l2Qb in the length-based expansion of Y n+ \ do not contribute to the 
first member of eq. (|27t for rj\/ = t + 1, (m'" = i) in accordance with eq. (fT2l . 
When instead m'" < i then m'" + 2 = (n + 1 — 2i) + 1 = m' + 1 since in such a 
condition m! = min(i, n + 1 — 2t) = n + 1 — 2t. Considering that when To = 
and sq = 1 the operator £\{n) may be represented as 



£y(n) = £l{n)£\(n) 



]J £ ( n - 2i J ]J £i ( n - 2r - j J 

io=o ^ io=o ^ ' 



(28) 



may be rewritten as follows: 



A(n)__ E z^rh 

(r,s)es , t r n _ 2t/Ti>0 



(f,S)€S; 



t,n— 2i/rj >0 



(29) 



where r runs from 1 to m'" + 1 in accordance with eq. (I20l ). We are interested in 
summing on r both members of eq. (1291 getting: 



E 



l nO r\ f>Ti rsi rT r ns r 
i L 'Q L -l L -Q L 'l ■ ■ ■ L '0 ^1 



m"'+l 

E 

r=l (r,S)e^ n _ at/T1>0 
m"'+2 

E E 

r=2 (r,5)65- n+1 _ 2VTi=0ai=1 



(30) 



As already discussed just before eq. (|28T ) if m'" = t then r may be stopped at 
m!" + 1 = m' + 1 whereas if mf" < t then m!" + 2 = m' + 1 as well. Thus we 
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are legitimated to write down that 



If] m"'+l 

£ £ A(n) E [£S 1 ^ 1 ---^£r]i 

t=0 r=l (t,s)6S; 



(31) 



,n— 2t/r 1 >0 

L-iH m'+l 

EE E [^r-.-^li- 

t=l r =2 (r,s)eST ... . . 

We have thus completed the demonstration of theorem 1 since we have shown that 
£o(n)Y n (£i(n)y n+ i) when Y n (Y n+ \) is expressed by eq. (fl2l generates all and 
only the ordered operators of Y n+ 2 as given by the same equation beginning with 
Q{n) withri > (n = 0). 

Remark 1. Exploiting a treatment analogous to the one used to demonstrate The- 
orem 1, it is possible to prove that the solution of eq. ([7]) given that Yq = A and 
Y\ = may be written down, for any n > 2, as follows 



rr «-2 i 

2 J 



Y n 



s r =0 



t=0 



A 



(32) 



where all the ordered contributions to be considered in the summation are those 
finishing with C\C\. 

To clarify the notation, we give the expression of I5 : 

The constraints extract from the first term only one ordered operator of length 2 
and from the second one, two ordered contributions only, of the same length 2: 

Y 5 = [A(3)A(2)A(1)A)(0) + £o(3)A(l)A(0) + £i(3)£ (2)A)(0)] A. 

Since the two Cauchy problems we solved give rise to two independent solutions 
of eq. £[]), its general solution, for any n > 2, may be written down exploiting eqs. 
C2]) and ([32]> 



Y n 



(AO L l }2/Tr = l,S r =0 

t=0 



rl"-i| 



^4+ E {4^1 

t=0 



t r n-l-2t 



}iB. (33) 



4 Applications 

In this section we show the effectiveness of eqs. (fT2l) and (l32l by exactly solv- 
ing a non-trivial example of eq. £[]) formulated in a complex vectorial space V 
of even dimension N. Following Dirac denote by \v) a generic vector of V and 
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by (v\w) := (w\v)* the scalar product between the two vectors \v) and be- 
longing to V. If an orthonormal basis B of V is composed of the N vectors: 
|1) ,|2), ...,[«},..., Ij), ... , | A^) such that = 5ij, let's introduce the tran- 
sition operator Ty := acting upon a generic vector \i') of B as follows: 
Tjj|i') = 6ji'\i). We define the following linear operators acting on V where the 
N/2 coefficients {ci,i = 1,2, ... , N/2} are complex numbers: 

N/2 N/2 

M+ = ^2aT iiN ^ +1 = Y^a\i)(N-i + l\ (34a) 
i=l i=l 

AT/2 

M_ = (A-/+) 1 " = ^ c*|iV - i + (34b) 
i=i 

N/2 

D+ = J2\c t \ 2 \i)(i\ (35a) 
i=l 

A 7 " 

D-= |cj\r-i+i| 2 |*)(*l ( 35b ) 

i=N/2+l 

M = D + - D_ (36a) 
D N = D + + D_. (36b) 
When \a\ = pi = p whatever i is it is easy to check the following properties: 

M±\v) = 0, y\v)eV (37a) 

[M+,M_]= M (37b) 
M ± Af T = D± (37c) 

M_D + = p 2 M_ => D + M + = p 2 M + (38a) 

M + D + \v) = => D+M-k) = 0, V|u) € V (38b) 

M_D_|u) = => D-M+\v) = 0, V|«) € V (38c) 

M + £>_ = p 2 M+ => D_M_ = p 2 M_. (38d) 

We point out that eqs. (I37ab - (l37cb and (138bl) - (138cl) hold without any restric- 
tion on the coefficients c n and that when N = 2 and the corresponding unique 
coefficient c± = 1 M + , M_ and Mq may be traced back to the well-known Pauli 
matrices. In the following analysis, we will make use only of the properties (I37ab - 
(I38db . rather than using the explicit representation given by eqs. (I34al) - (l36bb - It is 
therefore worthy to underline that generally speaking, other sets of M_ and D_ 
matrices fulfilling the properties (I34al i- (l38dl i exist. 
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The scope of this section is to exploit eqs. (fl2l and (l32l to provide the closed form 
of the following general Cauchy problem 

Y n+2 = C (n)Y n + £i(n)Y B+ i, n € N*, (39a) 

y =F , Yi=Fi, (39b) 

where Y n is a vector of V, Yq,Y\ are initial conditions fixed at will in V and 
the generally noncommuting coefficients operators Co(n) and C\{ri),n € N* are 
defined as functions of n, as follows: 

£o(n) = < ./ _ , , (40) 





71 


even 


M_ 


/? 


odd, 


M_ 


11 


even 




n 


odd. 



£i(n) = ^ ~ , , ' (41) 

[ M + n odd. 

To this end it is enough to consider the two initial conditions 

(I) Y = 0, Y\ = Y\. (42) 

(II) Y = Y , Y 1 =0 (43) 

4.1 Solving the Cauchy problem (I) 

Consider the operator [Cq 1 C^ 1 Cq 2 C S2 . . . Cq t C s {] \ and observe that as consequence 
of eqs. (I37ab . (I40l . (|4TI ) it vanishes as soon as one among the r integer exponents 
Ti is greater than 1. In this case eq. (fl2l indeed exhibits "£o- se g ment s"containing 
M£. Moreover if r > 1 then sm > 1 and the operator £f Q C\ 2 . . . £ft C\ r ]i 
vanishes because it includes the product 

. . . d(k - 2ri - sx + 1)£ (A: - 2r x - si) . . . (44) 

which coincides with M± in view of eqs. (l40l) and (|4TT) . Summing up, when the 
length r of the operator [Cq 1 C^ 1 Cq 2 £{ 2 . . . CQ r Cl r ]i exceeds 1, its ordered expres- 
sion given by eq. (HJ) vanishes. This property, of course, strictly related to the 
operators Cq, C\ greatly simplifies eq. (fT2l where the ordered values of t of inter- 
est become t = (for any n) and t = 1 (for n > 4). The reason is that t > 1 
and r = 1 requires n > 1 and then leads to a vanishing operator, as previously 
discussed. 

We thus concentrate ourselves on the evaluation of the expression of two oper- 
ators {C^C 1 ^ 1 }! and {£j£"~ 3 }i. For the first one we have r = 1 and 

{ClCT x }i = [jOPqC^ 1 ]! = U(n - 2) . . . A(0) = M* . . . M+M_ (45) 

where * is — (+) if n is even (odd). The second and third equalities stem from the 
application of eq. ® and eq. (|4TT) respectively. Taking into account eq. (I37cb and 
the constraints on the complex entries q, we finally get 

{r^r 1 }! = { C 2 ?, + n odd ' (46) 

L u 1 1 \ p n Z M_ n even. 
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When t = 1 (n > 4) , the maximal length tm compatible with t is tm = 2 and 



Then it is easy to check that 
[C\CT\ = £ (n - 2)d{n - 4) . . . A(0) = 



p n 4 -D+ n even, 



n n— 3 



M_ n odd. 



(47) 



(48) 



Substituting eqs. (l46l) and (|47"T ) into eq. (fl2l) . taking into account that all the con- 
tributions from any t > 1 vanish, the solution of the Cauchy problem (fT5T ) may be 
cast in the following closed and explicit form 



n = 0, 

Fi n = 1, 

M_yi n = 2, 

(p n ~ i D + + p n ~ 2 M^ )Y 1 n > 2, even, 

p n ^D + + M_ )Fi n > 2, odd. 



(49) 



In passing we note that, when pj = p= lforanyz,Y n = (Z) + +M_)yi, whatever 
n > 2 is which means that we obtain a veiy simple constant solution in this case. 



4.2 Solving the Cauchy problem (II) 

We make use of eq. (1321 observing that since for n > 2 

{4/:r 2 } 2/Tr=liSr=0 = [£r 2 4 = A(n - 2) . . . A(1)A)(0)] 2 = (50) 
In addition we have for any n > 4 

{£§£r 4 } 2/Tr=liSr=0 = [ClC^Clh = £o(n-2)A(n-4) . . . £i(l)£ (0) = 

(51) 

Thus the solution of the Cauchy problem (II) under scrutiny assumes the following 
form 

Y n = 0, 

L Y n = 2, 
n > 2. 



By direct substitution it is easy to check that this very simple sequence of vectors 
of V is the (unique) solution of the Cauchy problem (II). Exploiting eqs. (|49l and 
|), we write down the solution Y n of the general Cauchy problem as follows: 



Y, 



X? 

Yt 

Co^Yq + M.Yi 
(p n ~ 4 D + + p n - 2 M_)Yi 



n = 0, 
n = 1, 
n = 2, 
n > 2, even, 
n > 2, odd. 



(53) 
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Of course, interpreting the initial conditions Yq and Y\ as playing the role of 
"summation constants", eq. (1531) may be also refereed to as the general solution of 
difference equation (I39ab . 

In conclusion we stress that the target of this section that is to find the general 
solution of a selected difference equation exploiting eqs. (fT2l) (l32l and has been 
reached by eq. (|53T >. We do believe that solving this not trivial toy-difference 
equation represents a good test to appreciate the potential and practical value of 
the theoiy we have reported in the previous sections of the paper. 

5 Concluding remarks 

The applicative potentialities of eq. © goes beyond the even-wide context of 
operator difference equations. We indeed emphasize that the treatment does not 
exclude that the unknowns might depend on some continuous variables and that the 
operator coefficients £o(n) and C\{n) may act also upon the single components of 
Y n . To appreciate the value of this observation consider the following difference- 
differential equation: 

Y n+2 (t) = C (n)Y n (t) + A(n)Y n+1 (54) 

where the dot denotes the first time derivative of Y n (t), which in turn means to 
derive each component of Y n . The two generally noncommuting operators £q and 
L\ are supposed for simplicity only to linearly mix the component of Y n . Let's 
then introduce the following operator 

C Q {n) = C Q (n)^-; £i(n) = £i(n) (55) 

By definition Co(n) and d/dt commute and if Co(n) is represented by a matrix, 
then Co(n) is a formal matrix as well, whose entries (Co(n))ij are differential 
operators defined as follows: 

(£o(n))ij = [Co(n)]ij^ (56) 

where [Co(n)]ij are the Co(n) entries. With the help of this notation, eq. (1541 
formally coincides with eq. (Q]) whose resolutive formula given by eq. (T53T ) ap- 
plies also to difference-differential eq. (1541 . provided that the initial conditions are 
appropriately re-interpreted (each component of Yq, Y\ must be thought as func- 
tions of t). When £o(n) and C\{n) are defined in accordance with eqs. (l40l ) and 
(|4TT ) respectively, where the constant p = 1 is further adopted, the solution of the 
Cauchy problem I) (Yo(t) = 0, Y\ = Y(t)) can be cast as follows 

r M_Yi(f) n = 2, 
Y*{t) = I D + f t +M_ n > 2, even, (57) 
[ D + + M_f t n > 2, odd. 
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As a final remark, we wish to emphasize that the strategic and far-reaching value of 
the treatment of eq. CD leads to a resolutive formula achieved without requiring a 
priori the "mathematical nature"of the unknown "object" Y n . Since in the formal 
framework of a linear second order difference equations eq. CD incorporates further 
ingredients such as noncommutativity and is in addition compatible with a possible 
dependence of Y n on continuous variables, we believe that out treatment of eq. £[]) 
should attract the interest of researchers in many areas from physical mathematics 
and engineering to economy, biology, social sciences etc, where discrete and/or 
continuous variables-based mathematical modeling play a central investigational 
role. 
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